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ABSTRACT 

We use frequency-dependent position shifts of flat spectrum radio cores to estimate the kinetic power 
of AGN jets. We find a correlation between the derived jet powers and AGN narrow-line luminosity, 
consistent with the well-known relation for radio galaxies and steep spectrum quasars. This technique 
can be applied to intrinsically weak jets even at high rcdshift. 
Subject headings: galaxies: active — galaxies: jets — quasars: general 



1. INTRODUCTION 

Active Galactic Nuclei (AGN) are some of the brightest 
objects in the Universe. This makes them visible to huge 
distances, and therefore useful in applications ranging 
from studies of Galactic ionized gas (through interstellar 
scintillation: lLovell et al.l (|2008f) ). to measuring positions 
on Earth to centimetre precision (via the tec hnique of 
geode tic Very Long Baseline Interferometry; iMa et afl 
(|2009f )). to cosmolo gical applications (th rough their use 
as standard candles HWatson et al.l f|201 lh ) . 

AGN play a key role in galaxy formation and evolution. 
The evolution of black holes and galaxies are tightly cou- 
pled throughout cosm ic history (jMagorrian et alJ 11998b 
iHasinger et al.l 120051 ). and nuclear activity of black 
holes imparts significant feedback on the surround- 
ing gas through radiatively and mechanically driven 
outflows. Radio-emitting jets of relativistic parti- 
cles inflate cocoons of radio plasma, which in turn 
drive bow shocks through the host galaxy and be- 
yond (Kaiser fc Alexander! [l9 97V On large scales, these 
shocks are observed to uplift the hot gas present in 
galaxy clusters. Even once jet activity ceases, jet- 
inflated radio lobes can rise buoyantly t hrough clus- 
ter gas, transporting the gas outward (e.g. iFabian et al.l 
12003b iForman et al.ll2005h . 

AGN activity is intermittent, triggered by either radia- 
tively efficient accretion onto the black hole (the so-called 
"cold mode" ) ; or radiativ e ly ine fficient ( "hot mode" ) ac- 
cretion iBest fe Heckmanl (|2012f h Powerful radio AGN 
are typically associated with the radiatively efficient 
mode, and often exhibi t high-excitation ionization lines 
(jHardcastle et all 120071 ). This mo de of AGN trigger- 
ing is often associated with mergers (IShabala et al.l[20ll 
iRamos Almeida et al.l[20T2l ). and was likely prevalent in 
a younger, denser Universe. By contrast, the majority 
of low-redshift radio AGN are low power, and show no 
such strong emission lines; these are consistent with be- 
ing fuelled by stea dy cooling of hot X-ray emitting gas 
in galactic haloes (|Best et al.ll2005HPope et al.ll2012l ). 

From a theoretical viewpoint, AGN feedback is a cru- 
cial ingredient of all galaxy formation models, truncating 
(the otherwise) excessive star formation in massive galax- 
ies at late times, and ensuring that present-day ellipti- 



cals are "red and dead" (|Croton et al . 200G: B ower et al.1 
[2001. This feedback can come from radiative pres- 
sure (IFabian et al.ll200ll. quasar winds (jNesvadba et al.1 
[200l, je ts dShabala fe Alexanderll2009fl or a combination 
of these (jHopkins fc Elvisll2010l ). 

A number of authors have argued that AGN feedback 
signific antly affects the star f ormation histories of AGN 
hosts dKavirai et al.l 20111: lAntonuccio-Delogu fc Silk! 
|2TM iKauffmann et al.l 1 27)03) . iRawlings fc Jarvisl (|2004f ) 



suggested on energetic grounds that a single AGN may 
be abl e to expel gas f rom m ultiple galactic haloes. Re- 
cently, iShabala et al.1 (|2011| ) confirmed this observation- 
ally: shocks from the largest (hundreds of kpc-scale) ra- 
dio AGN were found to truncate star formation in dwarf 
galaxies that themselves may have never hosted an AGN. 
This potentially important new mode of feedback de- 
pends sensitively on the combination of AGN jet power 
and environment density. In order to form a large radio 
AGN, the radio jet needs to pierce the dense environment 
of the AGN host without being disrupted by magnetohy- 
drodynamic instabilities. Weak jets in moderate density 
environments, or moderate power jets in dense environ- 
ments, are easily disrupted within dense galactic cores 
(|Alexandenl200a ). and are thus unable to do large-scale 
feedback. Therefore, accurate measurement of AGN jet 
power is crucial to quantifying the effects of AGN on 
galaxy formation and evolution. 

Measuring kinetic AGN power is difficult. Dynamical 
models, utilizing observed sizes and lumino sities of radio 
AGN , can be used to est i mate jet powers (iKaiser et al.1 
[l99"7l IShabala et al.ll2008HAntognini et al.ll2012t) . These 
models rely on some knowledge of environment into 
which the radio lobes are expanding, and AGN age. A 
closely related cavity power method relies on measure- 
ments of the work done by the r adio lobes in inflating 
cavit i es in the X-ray emit ting gas (jRawlings fc Saunders! 
119911: iBirzan et all I2008D . This method also requires 
knowledge of both the density of the gas into which the 
lobes are expanding, and age of the radio lobes. Ages of 
synchrotron-emitti ng electrons can be estima ted from ei- 
ther radio spectra ([Alexander fc Leahvill987l ) or dynam- 
ical models, and t he two estimates typic ally differ by a 
factor of two (e.g. iMachals kTet al.l 120041 ). Furthermore, 
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the cavity technique is limited to nearby, low-power radio 
galaxies, and the derived jet powers may thus be unrep- 
resentative of the overall AGN population. Other sources 
of error in jet power estimates arise due to uncertainties 
in the exact gas density profile: for example, cavity mea- 
surements typically assume that the radio lobes expand 
into a constant density atmosphere throughout their life- 
time. 

In this paper, we apply a fundamentally different 
method for measuring AGN jet power to a sample of flat 
spectrum radio quasars. This method, first dev i sed b y 
iLobanovl (|1998f > and further discussed by lHirotanil (|2005f ). 
is based on the interpretation of the frequency depen- 
dent position shift of flat spectrum radio cores in terms 
of synchrotron self-absorption within the jet plasma. We 
compare the derived jet powers with narrow-line lumi- 
nosities. 

Throughout the paper, we adopt the concordance cos- 
mology of M = 0.27, fl A = 0.73, h = 0.71. 

2. MODEL AND ASSUMPTIONS 

This work is based on the standard model of self- 
absorbed radio cores, in which the unresolved, flat spec- 
trum "core" in VLBI images of radio quasars is un- 
derstood to be associated with the part of the jet at 
which the optical depth is unity. The core shift with 
frequency is interpreted as being due to the changing 
synchro tron self-absorption opacity in a freely expand- 



and q r = 1. Finally, the parameter qs = 2 
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cnro tron se n-absorp tion opacity m a ircely cxpand- 
jet (jMarscherll 19771 lKonigllll98lD . Using this model, 
a measurement of the core shift between multiple fre- 
quencies along with a number of reasonable assumptions 
allows an estimate to be made of the physical parameters 
of the jet plasma, and therefore, an estimate of the jet 
kinetic energy flux. Here we recount the relevant details 
of the model, a nd expand on the a nalysis presented in 
ILobanovl (|1998f ) and lHirotanH (|2005| ) in the context of jet 
p ower measurements. 

iSchwartz et al.l (|2006l Appendix B) provides a detailed 
derivation of the equation for kinetic energy flux in a rel- 
ativistic jet in terms of the physical parameters. The 
relevant quantities are defined as follows: (3 — v/c is the 
jet bulk velocity and T = (1 — the correspond- 

ing bulk Lorentz factor, A the jet cross-sectional area, 
n e the density of relativistic electrons/positrons, B the 
magnetic field strength, u e the energy density of rela- 
tivistic electrons/positrons and ub = B 2 /8ir is the mag- 
netic energy density in cgs units. We assume that the jet 
plasma may be a mixture of relativistically cold matter 
with pressure much less than rest mass energy density, 
and relativistic particles with pressure p e = | 
kinetic energy flux is then 
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The parameters q e and q r define the relative contri- 
butions from cold matter to the internal energy den- 
sity and rest-mass density, respectively. The parame- 
ter q e = |(1 + u p /u e ) where u p is the internal energy 
density of cold matter (i.e. protons). The parameter 
q r = (1 + Pp/ Pe) where p p and p e are the rest mass den- 
sity of cold matter, and relativistic electrons/positrons, 
respectively. In a purely electron/positron jet, q e = 4/3 



is the component of magnetic field perpendicular to 
the jet velocity, qs depends on the magnetic field ge- 
ometry, with < qs < 2. Throughout this paper we 
assume the magnetic field is dominated by the perpen- 
dicular component, so that qB ~ 2. 

2.1. Geometry 

We assume the quasar radio jet has a half opening 
angle <j>, and is viewed at an angle 8 (typically < 10°) 
to the line of sight. The corresponding Doppler factor 
is 5 = [r(l — /3 cos 6*)p\ The radial extent of the jet is 
described by r, while the transverse extent is given by 
R = ?*tan0 « r<p. 

2.2. Magnetic field and particle distribution 

Following iKonigll ()1981[ ) and ILobanovl (|1998| ) we 

parametrize the magnetic field as B(r) = _Bi(r/lpc)" a ; 
and the electron density as n e (r) = ni(r/lpc) _fc ; where 
B\ and n\ are the magnetic field strength and electron 
density 1 pc away from the central engine. Here, r is the 
radial distance from the central engine. Observations 
sugge st the core shif t scales as the inverse of f r equen cy 
{1/v; IKovalev et all (l200l : ISokolovskv etTaTl (l20ll 'l. 
This is indicative of conical jet structure and equipar- 
tition conditions, since parabolic jet structure produced 
due to external pressure gradients is likely to result in a 
significantly fla tter core position-frequency dependence 
(|Lobanovlll998f ). We therefore reasonably assume con- 
stant opening angle in our analysis. For such a jet, 
1 < a < 2, and the precise value depends on the magnetic 
field geometry: a — 2 for the component of magnetic field 
parallel to the jet, and a — 1 for the componen t of mag- 
netic field perpendicular to the jet. Following ILobanovl 
(fl998h we assume that the perpendicular component of 
the magnetic field is the dominant one, such that a = 1. 
Conservation of particles then implies 6 = 2. 

We assume a power-law energy distribution for the rel- 
ativistic electrons/positrons of the form N(^) = Ao7~ s . 
No represents the scaling of the particle energy distribu- 
tion, s = 1 — 2a and a is the optically thin synchrotron 
spectral index, typically in the range —0.5 < a < — 1. 



For a = 
depth as: 



2.3. Core shift 
1, b = 2, Hirotani (2005) gives the optical 
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Here, <5 is the Doppler factor; v is the observing fre- 
quency; a the spectral index; m, B\ are electron densi- 
ties and magnetic field strengths at 1 pc from the central 
engine; and all units are given in cgs except r which is 
in parsecs. C(a) is a tabulated function of the spectral 
index. Setting the optical depth r = 1 gives the location 
of the core as a function of frequency, r{v): 
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r(» = F(z, a, 5, 7min , 0, 6)nl La Bp- 1 i 
F = 1.759 x 10 7 (2.964 x 10 9 C(a " 



2tt(1 + z) 



(3) 



Equation [3] can be rearranged to give the difference 
between radii (in parsecs) at which the core is observed 
at two frequencies, 



(4) 



Equation |4] contains two unknowns, so in order to de- 
termine the physical parameters of the jet from a mea- 
sured core shift, we must make some assumption about 
the relationship between particle density (ni) and mag- 
netic field strength (Bi) at one parsec, which we discuss 
in the following section. 

2.4. Calculating the magnetic field strength under the 
assumption of equipartition 

iSokolovskv et al.l ()2011l ) find that most sources for 
which core shifts are easily measured have r(v) oc u , 
consistent with equipar tition of energy between particles 
and t he magnetic field (|Lobanovlll998t ISokolovskv et al.l 
120111) . We thus reasonably assume that the magnetic 
field is in energy equipartition with the radiating parti- 
cles (ub — u e ). Then, 
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where (7) is the mean Lorentz factor of the relativistic 
electrons/positrons, and for our assumed power-law dis- 
tribution 
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Hence, solving Equation [4] for B\ gives 
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2.5. Jet power 



Assuming a particular jet composition (electron- 
positron or electron-proton), spectral index a, and 
Lorentz factor limits 7 m j n and 7 ma x ; from Equation [1] 
we can write 



Q iet = KBfr'f3 
k = 7.2x 10 4 V 



g r r-i 1 
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Taking a = —0.5, Equation |7J] yields 
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where luminosity distance Dl is in Mpc, frequencies v\ 
and Vi in Hz, and core shift Ar mas in milli-arcseconds. 
The constants q e = 4/3, q r = 1 and qs = 2 for an 
electron-positron jet in which the magnetic field is per- 
pendicular to the jet, and Q(T,6, </>) = j-A 4> 3 ^ 2 is a 
term corresponding to jet orientation, opening angle and 
velocity. 

3. RESULTS AND DISCUSSION 

3.1. Correlation with narrow-line luminosity 

We use the above formalism to derive AGN jet pow- 
ers for a sample of flat spectrum quasars with measured 
core shifts in the literature, and compare the derived 
jet powers to their narrow-line l u minos ities. In a sem- 
inal work, iRawlings fc Saundersl (|1991|) found a corre- 
lation between narrow-line luminosity and jet power in 
a sample of 3CRR radio galaxies and steep spectrum 
quasars. For each source, these authors estimated jet 
power by calculating the work done on the surrounding 
gas by the radio lobes, and then dividing by the age of 
the radio source. This method assumes some coupling cf- 
ficie ncy between the jet an d the surrounding gas (factor 
k in IRawlings fc Saundersf ). Narrow- line emission results 
from the radiation field of the black hole accretion disk 
ionizing nearby gas clouds, and as such, a correlation 
between jet power and narrow-line luminosity is to be 
expected, since a higher accretion rate should result in 
both higher jet powers and stronger ionizing radiation 
field. 

From the literature, we selected quasars with measured 
core shifts, jet parameters and 0[II] and 0[III] line fluxes. 
Where multiple core shifts were available (as is the case 
for most of our sample), we calculated the frequency- 
independent core shift function = Ar mas ( „^ 2 V0 ) f° r 

each frequency pair, and adopted the median value. The 
jet power was then calculated using Equation [7J Narrow 
line luminosities were calculated via Lnlr = 3(3Lo[ii] + 
1.5L [in]) (jRawlings fc Saunders! 1 19911 ). 

Figure [T] shows the resultant jet power - narrow-line 
luminosity relat ion for an electron-p o sitron jet. Also 
plotted are the IRawlings fc Saunders! (|1991|) data con- 
verted to our cosmology. There is a strong correlation 
between the core shift-derived jet powers and narrow- 
line lumin osity, and this correlat ion is consistent with the 
results of [Rawlings fc Saundersl . Our slope of 1.0 ± 0.2 
agrees with the Rawlings & Saunders value of 0.9 ± 0.2. 
The scatter of 0.6 dex about the best-fit relation is also 
comparable with the Rawlings & Saunders value of 0.5 
dex. It is worth emphasizing that those authors used 
a completely different method (cavity jet powers) and 
a completely different sample (steep spectrum) of radio 
sources. 
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Figure 1. Core shift jet power - narrow-line luminosity re- 
lation (fi l led p oi nts) for a pair plasma. Core shifts ar e from : 
Lobanovl fi9 98). Kov alev et all l2u0l). ISokolovskv et al.l (pOlll l. 
O' Sulli van & Gabuzda (2009), Pushkarcv ct al. (200( g). Doppler 



factor s and beam i ng an gl es are from: [La htccmnaki fc Valtaoial 
l(T999l ). IXu et al.l lfl999Tl. IPushkarey et al. [20 03). Lin e lum i- 
nosities are f r om: Lawrence**et*lj!l i 19961 ). IWillott et al.l (119991 ) , 
Grhnes^taIJJ22Q4V Ope n black squares show the results of 
Rawlings & Saunders ( 1991) for a steep-spectrum sample of radio 
galaxies. The horizontal lines show conservative detection limits 
(0.1 mas between 2.3 and 8.4 GHz) in core shift at z = 0.5 and 
2 = 1; and the vertical lines are similar limits in narrow-line lumi- 
nosity. The detection limits at z — 0.1 are too faint to appear on 
the plot. Jet powers and narrow-line luminosities in the upper right 
quadrant of the plot should be observable throughout our sample 
volume. Lack of data in the top left and bottom right sections of 
the plot suggest strongly that the correlation is real. Note that 
the normalization of both our and Rawlings & Saunders relations 
depends on the assumed 7 ma x and 7 m i n cutoffs. 

3.2. Comparison with other jet power estimates 

A number of authors have used altern a tive t echniques 
to estimat e jet power. ICavagnolo et al.l ([2010D built on 
the work of lBirzan et alT i 20081 ) in reporting a correlation 
between cavity-derived AGN jet powers and radio lumi- 
nosities. [Willott et al. (1999) similarly use AGN radio 
luminosity as a proxy for jet power, based on a correla- 
tion between narrow-line and radio luminosities. These 
methods are useful, largely because radio luminosity is 
a relatively easy quantity to measure; and yield similar 
scatter (0.8 and 0.5 dex respectively) to our core shift 
technique. However, one must be careful in extrapolat- 
ing these results to large s amples. Dynamical m odels of 
radio sources, on which the IWillott et al.l (|1999D method 
is based, only apply to powerful Fanaroff-Riley (FR) 
type II radio sources. On the other hand, weaker FR 
I radi o sources domina t e AGN counts in the local Uni- 
verse (pest et all [20051 ; IShabala et all [20081) . Further- 
more, these models are sensitive to environment: a jet 
with a given kinetic power will appear brighter at ra- 
dio wavelengths if it propagates though denser gas. By 
contrast, the core shift method is independent of envi- 
ronment. 

The slope of the cavity je t power - integrated ra dio 
power relation reported by ICav agnolo ct al.l (120101 ) is 
0.75 ± 0.14 at 1.4 GHz, and 0.64 ± 0.09 at 327 MHz. 
At low frequencies the AGN radio luminosity will be 
dominated by lobe emission; while at higher frequencies 
flat-spectrum components such as jet core and hotspots 



Figure 2. Distribution of beaming parameter © for the complete 
sample of Pushkarev et al. (2009). Red circles represent quasars 
(51 objects), blue squares are BL Lacs and galaxies (16 objects). 
Shaded regions show Poisson uncertainties. 



will also contribute. This flat-spectrum contribution will 
be greater for powerful jets, resulting in a flattening 
of the Qjet~^radio slope at higher frequencies. These 
factors, along with environmental dependence, are also 
like ly to be responsible for the larger scatter (0.8 dex) in 
the ICavagnolo et al.l (|2010() relation compared with our 
method. Nevertheless, the slopes reported by these au- 
thors are broadly consistent with our results. 

3.3. Estimating core shift jet power 

iPushkarev et al.l ([20091 ) derived jet properties for a 
large sample (67) of flat spectrum ra dio AGN. This sam - 
ple is complete to 1.5 Jy at 15 GHz ([Lister et al.ll2009l) . 
and is thus representative of the flat spectrum radio AGN 
population. Using this data, we find that the distribu- 
tion of the function Q(T,8, <fi) is roughly log-normal for 
quasars (Figure [2]), which constitute the majority of the 
sample. Quasars have a mean of (logG) = —0.55 and 
standard deviation of cri og e = 0.37. Due to insufficient 
statistics (16 objects) it is difficult to say whether the 
distribution of for BL Lacs/galaxies is similarly log- 
normal. Values of (log 8) = —0.85 and cri og e = 0.42 are 
obtained under this assumption. 

Importantly, these distributions are narrow enough for 
core shift measurements to be useful predictors of jet 
power even if jet orientation parameters are not known: 
for a given narrow- line luminosity, 68 percent of the AGN 
should have jet powers that are within a factor of 2.5 
from the "best guess" values. This is illustrated in Fig- 
ure [3j where we plot the predicted jet power Qjet.approx, 
which is calculated using the mean logG values above, 
against narrow-line luminosity. Unlike Figure [TJ AGN 
with unknown values of T, 9 and <f> are also plotted. The 
strong correlation between these two quantities (slope 
0.8 ± 0.1, scatter 0.4 dex) suggests that core shifts can 
be used to estimate the kinetic energy flux of flat spec- 
trum AGN even if the jet velocity, viewing angle, and 
opening angle are not known. The best-fit relation for 
quasars, assuming an electron-positron jet in equiparti- 
tion for a = —0.5, is 
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Figure 3. Jet powers predicted using Equation \E\ Open symbols 
represent AGN with known jet parameters, plotted in Figure [T] 
Filled symbols are AGN for which beaming parameters are not 
known. Core shifts can be used to estimate jet power even if the 
jet parameters are unavailable. 
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where again luminosity distance Z?l is in Mpc, frequen- 
cies v\ and i>2 in Hz, and core shift Ar mas in milli- 
arcseconds. The median jet power is a factor of 2 lower 
for BL Lacs, for the same core shift and other proper- 
ties, due to the different distribution of jet parameters 
(as given by O; see Figure [2]). 

For an electron-proton jet, the rest mass term in Equa- 
tion [7]is important. Fo r normal cosmic plasma, q r rj 2200 
(jSchwartz et al.ll2006f ). and thus 
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The minimum Lorentz factor 7 m ; n in parsec-scale jets is 
highly uncertain, and whilst the minimum Lorentz factor 
in the hotspots of radio galaxies is up to 10 3 in some 
sources, 7 m i n in the jets prior to reaching the termi nation 
shock could be much lower (jGodfrev et al.1 12009). For 
7min = 10 and = 100 the estimated electron-proton 
jet power is a factor of 15 greater than the electron- 
positron case. This method could potentially be used to 
place constraints on the jet composition, and the value 
of 7min in parsec-scale jets. 

The estimation of jet powers from limited data, as 
given by Equations [5] and O is possible because the flat 
spectrum constraint necessarily limits the AGN sample 
to objects with a narrow range of viewing angle, and Y 
and </> have narrow distributions. 

4. CONCLUSIONS 

We have calculated AGN jet powers from multi- 
frequency core shifts in flat spectrum quasars. Wc 



find a strong correlation between jet power and AGN 
narrow-line luminosit y , con sistent with the results of 
iRawlings fc Saunders! (|1991l) for steep spectrum radio 
galaxies. The core shift method has strong predictive 
power, allowing accurate estimates of jet power to be 
made in the absence of detailed jet and AGN orientation 
parameters. 

The geometry of the Universe makes the core shift 
method a potentially very powerful one for measur- 
ing AGN jet powers, particularly in the era of new, 
high sensitivity instruments capable of very high astro- 
metric precision suc h as the Square Kilometre Array 
(jGodfrev et al.l l2012t) . The angular diameter distance 
turns over at redshift 1.6, making it in principle even 
easier to measure a core shift at redshift 3 than at red- 
shift 1.5. This is also true for other jet power measure- 
ment methods, such as dynamical and cavity models. 
However, the luminosity distance increases with redshift, 
making the detection of steep spectrum objects difficult. 
The advantage of the core shift method is in the fact 
that the flat spectrum nature of core-shifting quasars al- 
lows Doppler boosting to bring even intrinsically faint 
radio sources above the detection limit. Thus, the core 
shift method will be able to probe higher redshifts and 
lower luminosities than traditional jet power measure- 
ment methods. In addition to measuring AGN jet power, 
this technique is potentially useful in studying the assem- 
bly and evolution of black holes in the very early Uni- 
verse, and constraining the composition and the value of 
7min in parsec-scale jets. 
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